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We demonstrate a 'PT-symmetry-breaking chaos in optomechanical system (OMS), which features 
an ultralow driving threshold. In principle, this chaos will emerge once a driving laser is applied to 
the cavity mode and lasts for a period of time. The driving strength is inversely proportional to 
the starting time of chaos. This originally comes from the dynamical enhancement of nonlinearity 
by field localization in PT-symmetry-breaking phase (TTBP). Moreover, this chaos is switchable 
by tuning the system parameters so that a TT-symmetry phase transition occurs. This work may 
fundamentally broaden the regimes of cavity optomechanics and nonlinear optics. It offers the 
prospect of exploring ultralow-power-laser triggered chaos and its potential applications in secret 
communication. 
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Cavity optomechanics is a rapidly developing research 
field exploring the radiation-pressure interaction between 
the electromagnetic and mechanical systems [1]. Thanks 
to this nonlinear interaction, the optomechanical sys¬ 
tem (OMS) provides an alternative platform for im¬ 
plementing many interesting quantum [2-8] and classi¬ 
cal [9-17] nonlinearity phenomena. In particular, one 
can strong driving the cavity so that the OMS enters 
into a regime of self-induced oscillations [9-12], where 
the backaction-induced mechanical gain overcomes me¬ 
chanical loss. Further increasing the strength of the 
driving laser, chaotic motion emerges both in the op¬ 
tical and mechanical modes requiring no external feed¬ 
back or modulation [13-16]. It is useful for generating 
random numbers [18] and implementing secret informa¬ 
tion processing [19-21]. However, to apply chaos into the 
secret communication scheme requiring low-power opti¬ 
cal interconnects, the chaos threshold should be reduced 
dramatically [22, 23]. 

The notion of parity-time (TT) symmetry, initially 
proposed in quantum mechanics [24], attracts recently 
enormous attention in the field of optics [25-39]. It is 
known [24] that PT symmetric Hamiltonians can exhibit 
a real eigenvalue spectrum in spite of the fact that they 
can be non-Hermitian. These systems have the inter¬ 
esting property to undergo an abrupt phase-transition 
where the system loses its PT-symmetry. At the excep¬ 
tional point (EP) pairs of eigenvalues collide and become 
complex. Typically, the transition between PT symmet¬ 
ric phase (real spectrum) to spontaneously broken PT 
symmetry (complex spectrum) occurs as a parameter in 
the Hamiltonian (which somehow controls the degree of 
non-Hermiticity) is varied [25-27]. This phase transition 
has been demonstrated in synthetic waveguides and mi¬ 
crocavities, and can induce unique optical phenomena 
including loss-induced transparency [28], power oscilla¬ 
tions violating left-right symmetry [29], low-power opti¬ 
cal diodes [35] and single-mode laser [36, 37]. A natural 


question is whether it could influence the chaos dynam¬ 
ics (especially the chaos threshold) significantly. More¬ 
over, the crossover between the TT-symmetry theory 
and chaos in optomechanics remains largely unexplored, 
which may substantially advance the fields of cavity op¬ 
tomechanics and nonlinear optics. 

Here, we propose a TT-symmetry-breaking chaos [40] 
by investigating the nonlinear dynamics of a TT- 
symmetric OMS. In contrast to normal chaos, it is con¬ 
trollable by manipulating TT-symmetry phase transi¬ 
tion and allows an ultralow driving threshold, such as 
the weak-driving regime, i.e., the driving strength 
is smaller than the cavity decay rate 7. Physically, in 
TT-symmetry-breaking phase (TTBP) the field localiza¬ 
tion induces the dynamical accumulations of the optical 
and acoustical energy in the passive cavity, corresponding 
to an increasing optomechanical nonlinearity with time 
evolution. This ultimately leads to the result that the 
chaotic motion can be triggered even if an ultralow-power 
driving laser lasts for a period of time. As far as we 
know, this unconventional PT-breaking chaos, featuring 
an ultralow threshold, is identified for the first time in op¬ 
tomechanics, which may stimulate further investigations 
and applications in different optical, acoustic, or electric 
PT-breaking devices, with other types of nonlinearity in 
the loss element. 

System. — Consider an OMS coupled to an active cav¬ 
ity via optical tunneling, i.e., the TT-symmetric OMS 
[see Fig. 1(a)], with the Hamiltonian 

H = Hc +hujmb^b—hgQd\di{P-\-b), (1) 

where di (a|) and h (b^) are the annihilation (creation) 
operators of the passive cavity mode and mechanical 
mode, respectively. The passive cavity Si couples to an 
active cavity 0,2 (with the same resonance frequency Wc), 
and is driven with frequency uid, amplitude fid- Here 
fid is related to the input power P and decay rate 7 
by jOjjj = \J2P^/fuod. The Hamiltonian of cavity part 
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FIG. 1: (color online), (a) A schematic illnstration of a VT- 
symmetric OMS including a passive OMS (with cavity mode 
ai and a mechanical mode b) coupled to an active cavity 0,2 
with tunneling strength J. Here a)" and are the input 
and output of a driving field with frequency LJd', 7 ( 7 m) and k 
are respectively the optical (mechanical) decay rate and gain 
rate, (b) The realization of "PT-symmetric OMS in coupled 
microtoridal resonators, (c) The phase-diagram in terms of 
J /7 and n/'y indicates the VTSP and "PTBP in the passive- 
active system, as well as the parameter regime corresponding 
to passive-passive system. 

can be written as Hc/f^ = Ac(d|di - 1 - 0202 ) — J{a\a 2 + 
h.c.) — — oi), with Ac = Wc — in a frame ro¬ 

tating with uJd- The dimensionless position and momen¬ 
tum operators of the mechanical oscillator are defined 
as X = (l/-\/2)(6'^ -I- b) and p = (f/-\/2)(6'^ — b). The 
third term in Eq. (1) describes the radiation-pressure in¬ 
teraction between the passive cavity and the mechani¬ 
cal oscillator with coupling strength g [41]. In princi¬ 
ple, this T^T-symmetric OMS could be realized in many 
types of OMS [1], such as the coupled microtoroidal res¬ 
onators [35-37] [see Fig. 1(b)]. 

To explore the nonlinear dynamics of system, we em¬ 
ploy the semiclassical Langevin equations of motion (set¬ 
ting o = (d), o is an any optical or mechanical operator) 


X = ojmP, (2a) 

p = -i'jm/‘2)p-uJmX + V2go\ai\‘^, (2b) 

di = (—lAc — 7/2)01 -I- iJa 2 + iV2gQaiX -I- (2c) 

0,2 = {—i^c + K/2)a2 + iJai- (2d) 




FIG. 2: (color online). Power spectrum LnS(a;) of the intra¬ 
cavity field Ji = [flip in "PTSP (a) and "PTTBP (b). The inset 
of (b): bifurcation diagram versus J/ 7 . Lyapunov exponent 
versus (c) J/ 7 , inset of (c) K/ 7 , and (d) Hd /7 correspond¬ 
ing to a fixed time interval with 0.5 /rs. The shadowed blue 
region of the inset of (b) and (c) correspond to "PTBP, and 
the red arrows indicate the EPs. The inset of (d) corresponds 
to a case of weak driving, i.e., the black square in (d). Ac¬ 
cording to recent microcavity experiments [35], the parame¬ 
ters are lUc = 190 THz, tJm,/27r = 23 MHz, 7/271 = 1 MHz, 
go = 7.4 X 10 “® 7 , 7m/27r = O.O 387 , k = O. 87 , Ac = uj,n, and 
(a-c) P=lg,W (or Hd /7 ~ 4000), (d) J = 7 and J = O .27 
correspond to the blue circle and red square, respectively. 

T^T-symmetry phase transition when J /7 or k/j passes 
through the EP, i.e., J = (7 -I- k)/ 4, where the eigenval¬ 
ues and the corresponding eigenstates of system coalesce. 
This property might influence the nonlinear dynamics of 
system significantly. Here the EP is obtained by diag¬ 
onalizing the coefficient matrix of Eqs. (2c), (2d) under 
the condition of ignoring the optomechanical-interaction- 
induced shift of EP, which is valid in the weak-coupling 
regime 50/7 < 1 [42]. 

Governed by Eqs. (2), the evolution trajectory of sys¬ 
tem is specified by the initial condition. To charac¬ 
terize the stochastic properties of system, we linearize 
Eqs. (2) and introduce the evolution of a perturbation 

S = {Sx,dp,&iirjSaii,Sa 2 r,ba 2 i), i.e., S = M(5 with coeffi¬ 
cient matrix 


Here the quantum correlations of photon-phonon has 
been safely ignored in the semiclassical approximation, 
which is valid in the concerned weak-coupling regime, 
i.e., go/"f ^ 1 [13, 16]. Equations (2) show that the 
intracavity field intensity and the mechanical deforma¬ 
tion influence each other during evolution via the op¬ 
tomechanical interaction. Generally, this can lead to 
a chaotic motion of the optical and mechanical modes 
at a very high-power laser-driving [13, 16]. However, as 
shown in Fig. 1(c), the present coupled OMS features a 
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which characterizes the divergence of nearby trajectories 
in phase space. Here ba^r and ba^i are respectively the 
perturbation of the real part and imaginary part of Ok 
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FIG. 3 : (color online). Evolutions of (a) the intracavity in¬ 
tensities /i (passive cavity), I2 (active cavity), and (b) the di¬ 
mensionless mechanical displacement x in VTSP, i.e., J = 7. 
The insets: the trajectories of (a) the optical mode, i.e., the 
first derivation of /i versus 7 i and (b) the mechanical mode, 
i.e., p versus x, in phase space corresponding to a fixed time 
interval 8 /rs —> 9 ^s. The purple dashed lines indicates the 
amplitudes of Ii, I2 and x at 1.6 ps. The purple arrows indi¬ 
cate the approximate evolution direction of the trajectories. 
The parameters are same as Fig. 2 except for P « 0.02 pW, 
i.e., Q.d = 0.57. 


(A: = 1,2). 

VT-symmetry-hreaking-induced chaos. — By numeri¬ 
cally solving Eqs. (2), in Figs. 2(a) and 2(b) we respec¬ 
tively present the power spectrum of the field intensity 
1 1 = |aip in VT'SP and T^TBP, corresponding to a fixed 
driving power P — \ pW and time interval 8 —>■ 9/iS. On 
one hand, it shows that the cavity mode is periodically 
modulated by the mechanical oscillation with frequency 
uim in 'PTSP at a weak tunnelling coupling J = O.467 
[very close to the EP J = (7 -|- k)/4 = O.457]. When J is 
increased so strong that it could be comparable with (or 
larger than) the optical linewidth 7, a resolved normal- 
model-splitting with width Aw = j2 

appears around the center frequency and mechanical 
sideband due to the strong tunneling effect between the 
active and passive cavities. On the other hand, de¬ 
creasing J/7 pushes the system into the T^TBP, and 
the power spectrum becomes a continuum characteriz¬ 
ing the emergence of chaotic motion. Here the driving 
power 1 /rW is well below the threshold of normal op¬ 
tomechanical chaos for more than 3 order due to the held- 
localization-enhanced nonlinearity. Moreover, in the in¬ 
set of Fig. 2(b), we present the bifurcation diagram (cor¬ 
responding to the time interval 8 —)■ 9 ps) as the control 
parameter J/7 is varied, which shows that the system dy¬ 
namics experiences regular, period-doubling bifurcation 
to chaos behaviors from VTSP to T^TBP. 

In Figs. 2(c) and 2(d), we present the dependence of 
Lyapunov exponent [43] on J/7, /t/7 and ^d/l- The Lya¬ 
punov exponent is defined by the logarithmic slope of the 
perturbation Six [Six = [oi-I-5ai|^ — |ai|^) versus time t, 
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FIG. 4 : (color online). Evolutions of (a) Ii, I2, and (b) x in 
"PTBP, i.e., J = 0.27. The optical and mechanical trajectories 
in phase spaces are presented in the insets of (a,b). (c) The 
optical trajectory corresponding to time interval (hi). The 
inset of (c) presents the chaos dynamics in the long-time limit. 
The purple arrow in (c) indicates the approximate evolution 
direction of trajectory. The parameters are same as Fig. 3 . 


and characterizes the separation of trajectories for iden¬ 
tical systems with infinitesimally close initial condition. 
It is obtained by numerically solving Eq. (2) and the per¬ 
turbation equation <5 = M5 [13]. It clearly shows that the 
chaotic motion appears when one tunes the system pa¬ 
rameters so that VT symmetry is broken, which provides 
an alternative method to control chaotic dynamics with 
T^T-symmetry phase transition. Moreover, Fig. 2(d) also 
shows that the chaos allows an ultralow threshold, such 
as the weak-driving regime fid = O.Sy corresponding to 
P = 0.02 pW for microcavity system [35]. Note that, the 
oscillation amplitude of cavity intensity is not very large 
in the VTSP. However the semiclassical equation (2) is 
still valid, which is ensured by the bad-cavity regime, i.e., 
gxx/n ^ 1 and by using classical driving [13, 16]. Hence 
the Lyapunov exponent in Fig. 2 obtained from Eq. (2) 
can be trusted, which is also consistent with the evolu¬ 
tion of system Fig. 3. 

Dynamical trajectories in different VT phases .— To 
understand the ultralow threshold chaos in VTBP, in 
Figs. 3 and 4, we present the evolutions of the intracavity 
intensities Ji, 1 2 and the mechanical displacement x in 
VTSP and VTBP, respectively. 

Figures 3 show the amplitude oscillations in Ix, I 2 and 
X, characterizing the strong energy exchanges between 
the passive cavity and the active cavity as well as the 
mechanical oscillator (see the purple-dashed line). This 
ultimately leads to the periodic dynamics of the optical 
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and mechanical modes at a low driving-power in T^TSP. 
Accordingly, the optical and mechanical trajectories are 
limited into the regular circles with periodically varying 
radius with evolution (see the insets of Fig. 3). Physically, 
this is because, in T^TSP, the cavity-tunnelling effects 
characteristized by J is stronger than the intracavity lo¬ 
calization effects characteristized by (k -I- 7)/4. With¬ 
out strong optical driving, the weak intracavity intensity 
Ii cannot induce enough optomechanical nonlinearity to 
lead to a chaotic motion. 

Figures 4 show that, in T^TBP, the optical fields are 
localized respectively into their cavities (no clear energy 
oscillation between two cavities), and the chaotic mo¬ 
tions of the optical and mechanical modes appear af¬ 
ter a proper evolution even in the weak-driving regime, 
i.e., rid/T = 0-5. Comparing with the case of VTSP, 
here the evolution of Ii is more complicated, which can 
be divided into three distinguishable processes, (i) An 
analogous collapse and revival appears due to the inter¬ 
play between the optical loss and driving. The trajec¬ 
tories in phase space feature a set of limited circles ac¬ 
cordingly. (ii) It presents an exponential growth almost 
without oscillation with evolution, which demonstrates 
an enormous energy accumulation during this time in¬ 
terval. The approximate logarithmic-spiral-trajectory in 
phase space also indicate this property. Physically, it 
comes from the unidirectional energy transfer from the 
active cavity to the passive cavity in T^TBP, which also 
can submerge a probe laser injected from the active cav¬ 
ity and lead to a non-reciprocity transmission [35]. (iii) 
Chaotic motions of the optical and mechanical oscilla¬ 
tions emerges and lasts with evolution. The trajectories 
in phase spaces become very complicated [looks like a 
three-dimensional spirals, see Fig. 4(c)]. Accordingly, the 
mechanical trajectory goes trough similar processes [see 
Fig. 4(b)]. Note that the present chaos is bounded in the 
chaos regime due to the enhanced optomechanical nonlin¬ 
earity, as shown in the inset of Fig. 4(c), i.e., the evolution 
of system in the long-time limit. Based on the above, 
one can conclude that the essence of 'PT-symmetry- 
breaking chaos is the localization-induced dynamical- 
intensity-accumulation in PTBP. Therefore this chaos 
can be triggered even if an arbitrary small amount of 
power is applied to the system. Because strong intercav¬ 
ity intensity can be accumulated during the processes (i) 
and (ii) due to the effect of field localization. It can in¬ 
duce enough large optomechanical-nonlinearity required 
by chaotic motion. 

Comparison with the normal ehaotic dynamics. — The 
above PT-symmetric dynamics is quite different from 
the nonlinear dynamics in the normal passive-system, es¬ 
pecially in the chaotic regime. When k < 0 the PT- 
symmetric system returns to a passive-passive system 
[see Fig. 1(c)], whose dynamical evolutions are presented 
in Figs. 5(a) and (b). The comparison between Fig. 3(a) 
and Fig. 5(a) shows an enhanced optical tunneling phe- 



FIG. 5: (color online). The evolutions of Ii and I 2 in the 
passive-passive system when (a) P « 0.02 pW or Qd ~ O.Sy 
and (b) P ~ 500 mW or fid « 3 x lO^y. The insets present 
the corresponding optical trajectories in phase space. The 
black arrow indicates the approximate evolution direction of 
trajectory, (c) The approximate starting time r of chaotic mo¬ 
tion versus driving strength fld /7 in the active-passive system. 
The parameters are same as Fig. 2 except for (a,b) k = —O. 87 . 

nomenon in PT-symmetric system. Moreover, Fig. 5(b) 
shows that very strong optical driving (the order of mW 
for microcavity OMS) is required to obtain chaotic mo¬ 
tion in passive-passive system. The chaotic motion ap¬ 
pears at the beginning time and disappears with evo¬ 
lution due to the optical loss. Accordingly there is a 
transition from chaotic trajectory to regular trajectory 
in phase space, as shown in the insets of Fig. 5(b). Then, 
comparing Fig. 4(a) with Fig. 5(b), one can see that the 
PT-symmetry-breaking chaos is quite different from the 
normal optomechanical chaos. (1) It emerges after go¬ 
ing through two distinguishable processes, and lasts with 
evolution. (2) It allows an ultralow threshold (even is 
thresholdless), i.e., in principle, it can be triggered once 
a non-zero-power driving laser is employed and lasts for 
a period of time until enough nonlinearity is reached. 
The driving strength only resolves the starting time of 
chaotic motion, as shown in Fig. 5(c). Strong optical 
driving makes the chaos arise more quickly. 

Conclusions. — We have investigated the nonlinear dy¬ 
namics of an OMS coupled to an active cavity, which fea¬ 
tures a controllable PT-symmetry phase transition. We 
showed that, by tuning the tunnelling-to-loss or gain-to- 
loss rate, we can obtain an optomechanical chaos induced 
by the PT-symmetry breaking. Moreover we also pre¬ 
sented the dynamical trajectories of system in PTSP and 
PTBP, respectively. The increasing optomechanical non- 
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linearity with evolution in T^TBP makes the chaos fea¬ 
tures an ultralow threshold, which is required for many 
secret communication schemes. This study provides a 
promising route for controlling nonlinear dynamics, es¬ 
pecially the generation of chaos, with the concept of VT 
symmetry. It opens up new avenues for the study of 
ultralow-power-laser triggered chaos. Conversely, it also 
offers a method to distinguish the PTSP and PTBP via 
the experimentally detectable optical-spectrum. 
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